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Abstract. The medium-/ program of the Michelson Doppler Imager instrument on board SOHO provides 
continuous observations of oscillation modes of angular degree, /, from 0 to ~ 300. The data for the program 
are partly processed on board because only about 3% of MDI observations can be transmitted continuously 
to the ground. The on-board data processing, the main component of which is Gaussian-weighted binning, 
has been optimized to reduce the negative influence of spatial aliasing of the high-degree oscillation modes. 
The data processing is completed in a data analysis pipeline at the SOI Stanford Support Center to determine 
the mean multi plet frequencies and splitting coefficients. 

The initial results show that the noise in the medium-/ oscillation power spectrum is substantially lower 
than in ground-based measurements. This enables us to detect lower amplitude modes and, thus, to extend the 
range of measured mode frequencies. This is important for inferring the Sun’s internal structure and rotation. 
The MDI observations also reveal the asymmetry of oscillation spectral lines. The line asymmetries agree 
with the theory of mode excitation by acoustic sources localized in the upper convective boundary layer. The 
sound-speed profile inferred from the mean frequencies gives evidence for a sharp variation at the edge of 
the energy-generating core. The results also confirm the previous finding by the GONG (Gough et ah, 1996) 
that, in a thin layer just beneath the convection zone, helium appears to be less abundant than predicted by 
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theory. Inverting the nuiltiplet frequency splittings from MD!, we detect significant rotational shear in this 
thin layer. This layer is likely to be the place where the solar dynamo operates. In order to understand how the 
Sun works, it is extremely important to observe the evolution of this transition layer throughout the 1 1 -year 
activity cycle. 


1. Introduction 


Solar oscillations observed at the photosphere as fluctuations of Doppler velocity or intens- 
ity are usually described in terms of normal modes that are identified by three integers: 
angular degree /, angular order in, and radial order u. The frequencies of the normal 
modes depend on structures and flows in the regions where the modes propagate. These 
regions are different for different modes. In Figure 1 we show the spatial eigenfunction 
of a typical mode that propagates half way through the Sun in depth and at relatively 
low latitudes. The surface structure of the modes is represented by spherical harmonics, 
Y t m (0 , <f>). Because different modes sample different regions inside the Sun, by observing 
many modes one can, in principle, map the solar interior. 

The principal goal of helioseismology is to infer the structure and dynamics of the 
Sun’s interior from the oscillation frequencies of normal modes. The resolving power 
of these inferences depends on the accuracy of the frequency measurements, and on the 
frequency and angular degree bandwidths. Until now the frequencies of solar oscillations 
have been reliably measured only in a band around 3 mHz, from approximately 1.5 to 
4.5 mHz, where the mode amplitudes peak. These oscillations are acoustic (p) modes of 
low / and high n, or of high / and low n, and also fundamental (/) modes of high /. Modes 
of low / and low n have not been observed from the ground because their amplitude is very 
small, much less than 1 cm s 1 . Thus, one of the prime tasks of the MDI space experiment 
on board SOHO is to detect low-order modes. This requires long-term, continuous, stable 
observations free of atmospheric disturbances. 

MDI has three basic helioseismology programs: medium-/, low-/, and dynamics 
(Scherrer el al., 1995). The regions of the Sun covered by these programs are shown 
in Figure 2. The first two programs have been run with almost no interruptions since 
18 April, 1996. The medium-/ data are spatial averages of the full-disk Doppler velo- 
city out to 90% of the disk's radius measured each minute. This results in 26000 bins 
of approximately 10 arc sec resolution that provide sensitivity to solar p modes up to 
/ = 300. The low-/ observables are velocity and continuum intensity images summed into 
180 bins, with the intent of detecting oscillations up to / = 20. The dynamics program 
provides 1024 x 1024 images of the whole disk, thus covering all of the p modes up to 
l — 1500. However, the dynamics program can run continuously for only 2 months each 
year when the high-rate telemetry channel is available. Therefore, the frequency resolution 
of the dynamics data is somewhat limited. Also, the dynamics program cannot monitor the 
evolution of the solar structure and dynamics continuously. Therefore, it was important to 
extend the sensitivity of the medium-/ program to as high a degree as possible. 
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Figure I. The spatial structure of the solar mode of angular degree / = 20, angular order m — 16, and 
radial order n = 14. Red and blue show element displacements of opposite sign. The frequency of this mode 
determined from the MDI data is 2935. RK ± 0.02 // Hz. 
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Figure 2. Acoustic depth and radius of penetration of solar modes of frequency of 2.5. 3.5, and 4.5 mHz 
as a function of angular degree. / (black curves); and the areas covered by the MDI medium-/, low-/, and 
dynamics programs (shaded regions). The lower boundary of the convection /one is shown by the dashed 
line. The light-gray horizontal strip shows the location of the zone of the second ionization of helium. 
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The target for the medium-/ program is to observe the modes up to l = 300, because 
these modes have their inner reflection points in the He II ionization zone, and are, there- 
fore. particularly sensitive to the structure of this zone. Studying the He II ionization zone 
is important for measuring the helium abundance and calibrating the equation of state, 
and also for understanding the nature of supergranulation and the large-scale dynamics of 
the convection zone. 

To optimize the medium-/ observing program, we performed simulations of sever- 
al vector-weighted binning schemes, using a 20-hour series of full-disk Dopplergrams 
obtained on 25 January, 19%. An optimal set of Gaussian weights has been found that 
substantially reduces spatial aliasing in the angular degree range from 0 to 300 while 
preserving most of the power of the modes. 

In Section 2, we describe the technique used in the medium-/ program, the on-board 
data processing, and the basic parameters of the program. Analysis of the medium-/ data is 
presented in Section 3. In Section 4, we discuss initial inversion results for the sound-speed 
prolile and rotation. 


2. Medium-/ Program Technique 


The medium-/ data are transmitted through the low-rate telemetry channel (5 kilobits per 
second). The bit rate limits the spatial resolution of the observations. Prior to transmission 
the Doppler velocity is calculated from the original 1024x 1024-pixel CCD tiltergrams and 
binned on board the spacecraft. The binning can be performed in several ways depending 
on the scientific objectives of the program. For instance, Kosovichev (1992) considered 
spatially non-uniform binning schemes that maximize the range of the angular degree and 
maintain sensitivity to all the modes in this range. Alternatively, Gough (1992) suggested 
focusing on zonal and sectoral modes only, pushing the upper limit to I, = 1500. These 
and other possibilities will be explored in the course of the mission. However, the initial 
goal of the program is to obtain the most reliable measurements of modes of moderate / 
that provide information about most of the interior of the Sun. 

In observations performed with relatively low spatial resolution, the medium-/ data can 
be corrupted because of spatial aliasing. Spatial aliasing contaminates normal mode peaks 
in the power spectrum with the signal from high-degree modes. This poses significant 
problems for measuring parameters of the modal lines. Therefore, we have adopted the 
Gaussian vector binning scheme suggested by Milford and Schcrrcr ( 1992). This reduces 
the spatial aliasing by averaging the data with Gaussian weights centered on a regular grid 
(say. of 4 x 4- or 5 x 5-pixel spacing) on the CCD with a half- width approximately equal 
to the spacing of the grid. Such a 2-D Gaussian filter also reduces solar noise. We have 
evaluated Gaussian-weighted binning using a 20-hour long time series of one-minute full 
disk 1024 x 1024-pixel Doppler images obtained during the MDI commissioning and 
initial calibration observations. 
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l ii’tnv 3. Power spectra oT zonal modes computed using 20 hours of MDI full-resolution data box-car 
averaged onto a grid spacing of 4 x 4 pixels, covering the central part of the solar disk ir < 0.8/f) (/o/>) and 
for a 5 x 5 grid spacing covering 909f of the disk {bottom). 


2. 1 . Gaussian vixtor uinning 

For a vector binning scheme, the 2D filter function for a bin centered at the CCD pixel 
with coordinates (*«. jo) takes the form 


W (*», j n : i,j) = «(*«; 


( 1 ) 
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Tahle I 
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where / ) and are functions of the CCD column number, L and row number, 

/, respectively. While, in general, these functions can be arbitrary, the dyadic form of W 
is essential for the on-board processing. 

Following Milford and Scherrer (1992), we have represented u and v in the Gaussian 
form, i.e.. 

"('i>: i) = W o exp [-(< - *o) 2 /« 2 ] , (2) 


and similarly for /’(/o; /). Flere, a is the characteristic width of the Gaussians, and H o is a 
scaling factor. The centers (/ 0 , /<)) of the bins were uniformly distributed on a 5 x 5-pixel 
rectangular grid in the area within 0.98 solar radii of disk center. The total number of bins 
was 29804. Results were evaluated for a=2.5, 3.5, 4, 4.5, 5.4 and 6. For comparison we 
also computed oscillation power spectra for 5 x 5-pixel box-car averages over the full-disk 
and for 4 x 4-pixel box-car averages in the central area restricted to r < ().8/£. Each case 
matches the telemetry bandwidth limit. 

The data binning reduces the shot noise that determines the precision required in the 
on-board accumulations. We have estimated the coefficient of reduction of the shot noise, 
relative to the single pixel value as 


" "Ei.jH'li../) ■ 

The coefficients for the various values of a are shown in Table 1. 


(3) 


2.2. Power spectra of the spatially binned data 


Power spectra of the binned data were computed using two different techniques: the stand- 
ard method of remapping the images onto a rectangular grid of heliographic coordinates, 
and a method of direct projection on eigenmode masks (Kosovichev, 1992). The results 
obtained by both techniques are consistent with each other. 

The power spectra of zonal modes (m = 0) computed with 4x4 and 5 x 5 box-car 
averaging are shown in Figure 3. Obviously, the data sampled on the finer spatial grid 
resolve p modes at higher / better. Also, the relative amplitudes of the aliased modes 
crossing the original mode ridges are substantially smaller. However, observations with 
4 x 4-pixel averaging can only cover a smaller portion of the solar disk, increasing the 
spatial leakage from adjacent modes. This makes the mode ridges broader and complic- 
ates identification of individual oscillation mode frequencies. Therefore, for the initial 
observations, we require a 5 x 5-pixel grid which covers 90% of the full disk. 
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Figure 4 shows the power spectra of the sectoral modes (/// = /) computed with 
Gaussian weighting with widths a- 2.5. 4, and 5.4. As expected, increasing a reduces the 
amplitudes of the aliased ridges, hut it also reduces the amplitude of the original modes at 
higher/. Fora > 4, the amplitude of the ‘true' modes in the range / = 200 - 300 becomes 
too small. It seems that a — 4 represents a reasonable tradeoff between anti-aliasing and 
high-/ cut-off. In this case, the aliasing is substantially reduced up to / = 300. 

Therefore, for the initial medium-/ program observations, we have chosen to use the 
Gaussian vector-weighted data binning on a rectangular grid of 5 x 5 pixels, in the form 
of Hquations ( 1 ) and (2) with a — 4. The shot noise of the binned data is reduced to 0. 1 of 
the single-pixel value (Table I). If the noise per pixel per 60 s is 12 m s 1 then the noise 
per bin is 1 .2 m s 1 . The binned images are generated on board once per minute. 


3. Analysis of the Medium-/ Data 


3.1. POWHR SPHCTRA 

A standard helioseismic data analysis procedure (e.g., Hill et al ., 1996) has been applied 
to obtain the oscillation power spectra from the medium-/ data. 60 days of data have been 
processed. The images were apodized with a cosine bell between 0.85// and 0.89// and 
remapped onto a grid with 768 points around the equator and 128 intervals between the 
equator and pole. After this, they were Fourier transformed in longitude and multiplied 
by the associated Legendre functions /y"(ft) in colatitude ft to do the spherical harmonic 
decomposition. 

The power spectra of the sectoral and zonal modes, and also the ///-averaged power 
spectrum, obtained from 10 days of the data are shown in Figure 5. The ridges in the 
diagrams correspond to modes of different radial order n. The lowest weak ridge is the / 
mode. At low /, the ridges are not resolved in the plot of the / — // diagrams because the 
spectrum becomes too dense. 

The aliasing effect is strongest for zonal (m — 0) and sectoral (m = /) modes because 
these modes have the smallest spatial scale. It is clearly seen in the power spectra of these 
modes for / > 220. In the ///-averaged spectrum (Figure 5, bottom), the effect is much 
weaker. It is important that aliasing of the fundamental modes, /, and first-order acoustic 
modes, p\. be negligible up to / = 300. This makes it possible to measure frequencies of 
individual modes along these two lowest ridges in the / - v diagram. For higher p-mode 
overtones, individual modes are not resolved at / > 220 because the frequency difference 
between the modes of adjacent / becomes smaller than the line widths of the modes. In 
this case, instead of measuring individual peaks, we plan to use ‘ridge-fitting' techniques. 

The medium-/ program also provides good sensitivity to low-degree modes (Figure 6), 
thus covering almost the w hole range of the solar oscillation spectrum, in which individual 
p modes can be resolved. 
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Figure 5. Power spectra (/ — o diagrams) obtained from 10 days of the MDI medium-/ data for the sectoral 
(/ = w) modes {top), zonal (m = 0) modes {middle), and for the modes averaged over m {bottom). 
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(>. Power spectrum of/ = 0, I , and 2 modes obtained from 10 days of the MDI data. 


3.2. RESOLUTION OF low-frequency modes 

An important difference between the MDI medium-/ program and the ground-based 
networks is the ability to detect the low -n low-frequency modes that carry substantial 
independent information about the solar structure and dynamics. The lifetimes of these 
modes are much longer than the lifetimes of high-frequency modes. Therefore, in prin- 
ciple, the frequencies of these low-//, modes can be determined more accurately than the 
frequencies of the other modes. However, the amplitude of the low-// modes is very small, 
so long stable lime series are required to detect them. Two examples of low-frequency 
modes: the p mode of n = 1 and / = 50, and the / mode of / = 95, are shown in Figures 
7 and 8. Parameters of these modes have been measured for the first time. With longer 
time series, we should be able to observe more low-frequency modes, with even lower 
amplitudes. So far, the smallest mode amplitude we have been able to detect was about 
1 mm s 1 . 

3.3. Line asymmetry 

Because of the significant reduction in the noise, the medium-/ data have revealed inter- 
esting characteristics of the line profiles of the oscillation power spectra. Figure 9 shows 
the power spectra of the modes for / = 100, 200, and 300. The most interesting feature is 
the asymmetry of the line profiles. Though the asymmetry has been noticed in the ground- 
based data (Duvall et <//., 1993), frequencies of solar modes are usually determined by 
assuming that the line profile is symmetric and can be fitted by a Lorentzian, which would 
be the case if the solar// modes were damped simple harmonic oscillators excited by a 
stochastic source. However, this leads to systematic errors in the determination of fre- 
quencies (Hill et ///., 1996; Abrams and Kumar, 1996). Several authors have studied this 
problem theoretically and have found that there is an inherent asymmetry whenever the 
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Figure 7. Power spectrum of jpi mode of / = 50 obtained from 2 months ol the MDI medium-/ data. Modes 
of adjacent / leak into the spectrum because of imperfect spatial tillering. 



Figure 8. Power spectrum of / mode of / = 95 obtained from 2 months of the MDI data. 


waves are excited by a localized source (Gabriel, 1992, 1993, and 1995; Kumar etui 1994; 
Roxbourgh and Vorontsov, 1995). Physically, the asymmetry is an effect of interference 
between an outward direct wave from the source and a corresponding inward wave that 
passes through the region of wave propagation (Duvall et at 1993). Figure 10 shows a 
theoretical power spectrum of p modes of l = 200 obtained by Nigam and Kosovichev 
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Figure 10. Theoretical power spectrum of / = 200 mode for a model w ith a source of the uniform spectral 
density localized 80 km beneath the photosphere. The background noise is assumed to be uniform with the 
power spectral density 10“ 2 nrs ~ Hz’. 


(1996). This model is in good qualitative agreement with the observations. The degree 
of the asymmetry depends on the relative locations of the acoustic sources and the upper 
reflection layer of the modes. This opens the prospect of using the observations of the line 
profiles of solar modes to test theories of excitation of solar and stellar oscillations and of 
their interaction with turbulent convection. 


3.4. Frequency measurements 


The frequencies of the solar oscillations have been measured as described by Sehou ( 1992). 
The frequency splitting of mode multiplets due to rotation and asphericity is approximated 
using a polynomial expansion similar to that of Duvall, Harvey and Pomerantz (1986): 

6 

Vfilm ~ u nl + X] (4) 

k~ 1 

where u n ( is the mean frequency of a mode multiples and are orthogonal poly- 

nomials of degree k defined by 

i 

VP(1) = L and vP(m)pf ] (m) = 0 for i ^ j. (5) 

711 =-/ 


The polynomials, can be expressed in terms of the Clebsch-Gordan coefficients, 

syim 

L k(Unr 


r[ l] (m) = 


ww-km+k + ni , m 

{2l)W2l + I m "' 


( 6 ) 
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f it>nre II. Mean frequencies of the mode multiplcts ohtained from 2 months of the medium-/ data. The error 
bars show the standard errors multiplied by 1000. 


Since 

i 

Y P^\m) = 0 for A: > 0. (7) 

Hi - / 

the mean frequencies, /v„/, depend only on the spherically-symmetric component of the 
solar structure. Frequency splitting coefficients, for even A: depend on the aspherieal 
component of the structure, while the coefficients for odd A; measure the rotation rate (e.g., 
Gough, 1993). 

For the initial frequency measurements, we approximated the line profiles by Lorent- 
zians and used a maximum likelihood method to determine the parameters of the Lorent- 
zians, taking into account possible overlaps of lines. The effect of the line asymmetry 
was not taken into account in these initial measurements. The mean frequencies of mode 
multiplcts are shown in Figure I 1. The error bars in this figure indicate estimated errors 
multiplied by 1000. 


4. Initial Inversion Results 

4. 1 . Radial stratification 

We have determined the spherically-symmetric structure of the Sun by using the optimally 
localized averaging techniques (e.g., Gough, 1985; Pijpers and Thompson, 1992) to invert 
the mean frequencies of split mode multiplcts, v n \. Figure 1 2 shows the relative difference 
between the square of the sound speed in the Sun and a reference solar model. The 
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Figure 12. Relative differences between the squared sound speed in the Sun and a standard solar model as 
inferred from 2 months of the MDI data. The horizontal bars show the spatial resolution, and the vertical bars 
are error estimates. 


reference model was model S of Christensen-Dalsgaard et ai (1996). This model was 
a standard evolutionary model computed using the most recent information on nuclear 
reaction rates (Bahcall and Pinsonneault, 1995) and radiative opacity and the equation of 
state (Rogers and Iglesias. 1996). The gravitational settling and diffusion of helium and 
heavier elements were taken into account following the theory by Michaud and Proffitt 
(1993). 

The inversion results show that the maximum difference in the square of the sound 
speed between the model and the Sun is only 0.4%. Nevertheless, this difference is very 
important for understanding solar evolution and physical processes inside the Sun. Two 
features of the sound-speed profile are particularly notable. The first is a narrow peak 
centered at 0.67 R, just beneath the convection /one. This peak was previously detected 
in the LOWL (Basu et ai , 1996) and GONG data (Gough et ai, 1996) and is most 
likely due to a deficit of helium in this narrow region. The deficit of helium decreases the 
mean molecular weight and thus increases the sound speed. The deficit of helium could 
result from additional mixing of the material in the layer with the surrounding plasma if 
turbulence is generated in this layer because of rotational shear. Indeed, as we show in the 
next section, there is a strong radial gradient of the rotation rate in this layer. 

Another interesting feature is the sharp decrease of the sound speed compared to 
the model at the boundary of the energy-generating core, at 0.25 R. This variation is 
not explained by the theory. Kosovichev and Fedorova (1991) noticed a possible sharp 
increase of density at the edge of the core when inverting the BBSO data (Lihbrecht, 
Woodard, and Kaufman, 1990) if the spatial resolution of the inversion was increased at 
the expense of substantially increased errors of the estimates. The evidence from the MDI 
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data is more convincing. It is quite possible that the drop in the sound speed results from 
an overabundance of helium at the edge of the solar core. 

The steep increase of the sound speed towards the solar center can be explained if 
helium is less abundant than in the standard solar model. This indicates that the helium 
abundance profile seems to be more flat in the solar core than it is in the standard solar 
model. Gough etal. (1996) came to a similar conclusion from the GONG data. In principle, 
this variation of the helium abundance could arise from errors in the rates of the nuclear 
reactions or element diffusion. It could also be explained by uncertainties in the equation 
of state and in radiative opacity. Helioseismology is unable to distinguish among the 
possibilities without additional assumptions about the physical processes. However, if 
the transition at the edge of the core is really as sharp as we have found from the initial 
MDI data, then it strongly suggests material redistribution by macroscopic motions in the 
core, possibly induced by the instability of 3 He burning, as first suggested by Dilke and 
Gough (1972). Detailed studies of this phenomenon are very important for understanding 
stellar evolution. The deviations from the standard model of the spherically averaged 
structure of the solar core are fairly small; too small to solve the solar neutrino problem 
directly. However, if an instability really occurred, as predicted by theory (Christensen- 
Dalsgaard, Dilke, and Gough, 1974; Boury et al ., 1975; Shibahashi, Osaki and Unno, 
1975; Kosovichev and Severny, 1985; Merryfield, Too m re and Gough, 1991), then the 
basic assumptions of solar evolution theory about the thermal balance in the core must be 
questioned. This is certainly important for understanding the nature of the solar neutrino 
problem and for using the Sun as a laboratory for particle physics. 

4.2. Internal rotation 

The internal rotation rate at three latitudes, 0 (equator), 30 and 60 deg, is shown in 
Figure 13. The rotation is inferred from the odd coefficients of Equation (3) using 
a regularized least-squares technique (e.g., Tikhonov and Arsenin, 1977). The inversion 
results confirm the previous findings that latitudinal differential rotation occurs only in 
the convection zone, that the radiative interior rotates almost rigidly, and that there is 
a thin shear layer near the surface (Goode et ul . , 1991; Thompson et al . , 1996). The 
most interesting feature is that the transition layer (tachocline) between the radiative and 
convection zone is mostly located in the radiative zone, at least at the equator where it is 
also fairly thin, certainly less than 0.1 R. The layer seems to be wider at high latitudes. 
Previously, the parameters of the tachocline were estimated by Kosovichev (1996) and 
by Charbonneau et al. (1996) by fitting an analytical model to the data from the BBSO 
(Woodard and Libbrecht, 1993), and from the LOWL instrument (Tomczyk et al ., 1995). 

The sharp radial gradient of the angular velocity beneath the convection zone strongly 
suggests that the sharp narrow peak of the sound speed at 0.67 R seen in the structure 
inversion (Figure 12) is indeed due to rotationally-induced turbulent mixing in the tacho- 
cline (Spiegel and Zahn, 1992). The location and the width of the tachocline are also 
consistent with the requirement of recent dynamo theories of the solar cycle (Rudiger and 
Brandenburg, 1 995; Weiss, 1994). It will be extremely important to observe the variations 
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r/R 

Figure Li. Solar rotation rate inferred from 2 months of MDI medium-/ data as a function of radius at three 
latitudes, 0°, 30°, and 60°. The formal errors are indicated by the shaded regions. The arrows indicate the 
Doppler rotation rale directly measured on the surface. 


of the angular velocity in the tachoeline through the solar cycle in order to understand the 
mechanism of solar activity. 


5. Conclusions and Perspectives 

The MDI medium-/ program is designed for continuous monitoring of the structure and 
dynamics of the Sun’s interior. The initial results show that the noise in the medium-/ 
oscillation power spectrum is substantially lower than in the ground-based measurements. 
This enables us to resolve lower amplitude modes and, thus, to extend the range of 
measured mode frequencies. This is important for inferring the Sun’s internal structure 
and rotation. The MDI observations clearly show the asymmetry of oscillation line profiles 
in agreement with the theory of mode excitation by acoustic sources localized in the upper 
convective boundary layer. The sound-speed profile inferred from the mean frequencies 
gives evidence for a sharp variation at the edge of the energy-generating core, which is 
not accounted for by the standard evolution theory. The results also confirm the previous 
finding from GONG data that helium appears to be less abundant than predicted by 
theory in a thin layer just beneath the convection zone (Gough et ai , 1996). Inverting 
the multiplet frequency splittings from MDI, we have found significant rotational shear 
in this thin layer. This shear flow probably generates turbulence that mixes the plasma in 
the upper radiative zone. This layer (‘tachoeline’) is likely to be the place where the solar 
dynamo operates. 
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The first results of the medium-/ program suggest the initial binning scheme should be 
maintained for the first year of SOHO observations to obtain the time series sufficient to 
resolve low-frequency p modes. After the first year, it may be interesting to explore other 
possibilities, such as decreasing the remaining aliasing for modes of / = 200-300 to allow 
identification of a greater number of modes, or using an inhomogeneous binning scheme 
to focus on some particular groups of modes. Certainly, the decision should be taken 
after consideration of the performance of the other SOHO and ground-based helioseismie 
programs. However, the principal goal of the medium-/ program will remain the same - 
continuous measurement of the solar structure and dynamics as they change through the 
solar cycle. This is crucial for understanding the mechanisms for solar variability. 
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